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SUMMAET 



Proceeding from the basic assumptions of the Batho- 
Bredt theory on twisting failure of thin-walled columns, 
the discrepancies most frequently encountered are analyzed. 
A generalized approximate method is suggested for the de- 
termination of the disturbances in the stress condition of 
the column, induced by the constrained wg i nl e ltn g in one of 
tho ond sections. •"••'/••"J 

INTRODUCTION 



The evaluation of the stresses induced by the appli- 
cation of a twisting moment in a structure consisting es- 
sentially of thin-walled cylindrical tubing and of arbi- 
trary straight section, is generally carried out with the 
simple and easily applied formulas given in tho well-known 
Batho-Brodt approximate theory. 

The basic assumptions of this theory are as follows: 

la The wall thickness of the tube is thin enough to 
render the variations in unit stress along the 
perpendiculars to the walls negligible. 

2. The shape of the straight section of the tube is 

preserved during the strain induced by the ap- 
plication of the twisting moment. 

3. The sections of the tube Including those at the 

' end are bo restrained as to offer no resist- 
ance to w r Vn^llng f. the displacements of points 
of the. tube in direction of the generating axes 
of the cylindrical surface boing freely per- 
mitted. 



* "Sulla. Torsione di Tubi a Paroti Sottili Perf ettamente In- 
castrati ad una Estremita." L 'Aerotecnica, vol. XVII, 
no. 10, October 1937, pp.. 850-861. 
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TChilo the first tiro assumptions may he held to he 
confirmed with sufficient approximation in the majority 
of practical cases, the third is far from it - at least 
one of the end sections heing joined to the other much 
moro rigid structure or otherwise restrained in such a way 
that its wrinkling is not permitted. The constrained 
wrinkling in one of the end sections of the column ohvi- 
ously produces stresses in tho diroction of the generat- 
ing axos of tho column, the values of which certainly 
cancol in tho froe-ond section. 

On structures whose form would introduce considcrahlo 
wrinkling, such as those with considerahly elongated sec- 
tions, the previously cited stresses may assume dangerous 
proportions in certain points, especially if accompanied 
hy stresses due to other concomitant causes ("bending, for 
example ) • 

Hence the importance, in certain cases, attaching to 
tho evaluation - oven approximate - of the stresses in- 
duced hy constrained wrinkling in one of the end sections 
of the structure. For convenience of treatment, -i-t is ad- 
mitted that the stress conditions created when applying a 
torque on a column with one end perfectly restrained, may 
he considered as the sum of the stress condition as de- 
fined hy the Batho-Bredt theory and a secondary stress 
condition, which nullify in the restrained section the 
wrinkling pertinent to the primary stress condition. 

Tho secondary unit stresses at any section of the 
column must have zero resultant and zero resultant moment. 

According to the Batho-Bredt approximate theory, uni- 
tary parallel tengontial stresses exist in the plane of a 
gonorio section of the column, normal to the generators, 
in any part of the median contour of the section, and hav- 
ing the value: 




where M is the applied torque, S the area enclosed hy 
the median contour of the section, and 6 the wall thick- 
ness at the point in question. 

It should he noted that the section shape has no ef- 
fect on the value of T. This, however, does not hold for 
the displacements of the points of the column Induced hy 
that stress condition. 
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From the Batho-Bredt theory. It is also seen that all 
straight sections of the. column- wrinkle- equally, and turn, 
with tho rotation of a rigid body, in relation to each 
othor . 

Tho coordinates of a point of the column wall are: 

The distance z from the plane of one of the end 
sections which we suppose to he perpendicular to the gen-* 
orators of the walls; 

The length a of the arc of the median contour of 
the straight section to which the point in question he- 
longs, contained between it and any generator assumed as 
origin. 

To simplify the writing of the formulas, the unit 
length is chooen so that the perimeter of the median con- 
tour of the section of the column is equal to 2tt. 

Let w and t represent the components of tho dis*« 
placement of a generic point of tho wall with respect to 
z and with respect to the tangent to the median contour 
of tho straight section of the column at the point under 
consideration, visualized in the position preceding the 
application of torque. 

Tho relation of gonoral character is verified in each 
point: 

|Z+|*.I (2) 
d s a z G- 



whore 0 is the modulus of tangential elasticity of tho 
mat'orial. Let 0. he tho angle through which a generic 
soctlon has rotated in relation to its initial position 
undor tho effoct of the 'Batho-Bredt stress condition; fur- 
ther, lot r represent tho length of the segment con- 
tained botweon a generic point of the wall and tho point 
0 about which tho straight soctlon of the column, to which 
tho considered point bolongs, has rotated. Lastly, lot a 
roprasont tho anglo formed by this segment with tho tangent 
to' tho wall at tho point in question;" 

Prom what has already boon said rogarding tho naturo 
of the strain inducod by tho stress condition, according 
to "the Bntho-Bredt theory, follows (fig. l): 
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at 

dz 



r sin a ^ 



(3) 



where 



4G S 




(4) 



o 



as the cited theory stipulates. 

By establishing in the plane of the straight section 
of the column a system of orthogonal cartesian axes x; 
y, and z, and assuming as positive direction on the arc 
s, that which "brings axis x on axis y, equation (3) 
may he written in the form: 



whore a, b are the coordinates of 0 with respect to 
the chosen reference. 

The coordinates a, b cannot he determined on the 
basis of the Batho-Bredt theory, inaemuch as this, as has 
"beon shown, prescinds from the shape of the straight sec- 
tion of the column, upon which a and b depend. 

Tho position of point 0, about which the section 
in fact rotates, depends in roality on the limiting con- 
ditions which must comply with the condition of strain. 

Those requirements of the Batho-Bredt approximate 
theory are insufficient to determine the position of 0. 
For different positions of 0, then, we can always ob- 
tain the 3tress condition defined by that theory. From 
among those, any one may be chosen, imposing on the strain 
condition wholly arbitrary supplementary conditions which, 
however, oust conform to the assumptions upon which the 
above-nentioned theory is based. Thereupon also, essen- 
tially are based the various attempts made by several au- 
thors to determine the center of twist of columns in tor- 
sion patterned after the Batho-Bredt method. 

Here the supplementary conditions which are to satisfy 
the principal condition of strain are prescribed by the 
general consideration of the method adopted for the rep- 



ft = M (x _ a) |z _ ( _ „) 5£l 
dz dz L ds dsj 



(3a) 
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re sen tat ion of the total stress condition; a and h, how- 
ever, assume values which nullify- the components of the 
resultant handing moment of the secondary stresses acting 
along the column generators which reflect the end section, 
twisted under the effect of the Batho-Bredt stress condi- 
tion, in the plane in which it is restrained. 

Substituting the expression supplied "by equation (3a) 
in equation (2), we have: 



dw 
ds 



I + £! [" (y-b) ^ - ( X -a) fZ] 
G ds L dB dsj 



(5) 



whore 



dw 
ds 



ia givon in function of the Known quantity sup- 



plied by oquation (l), equation (4), and by the conforma- 
tion of tho straight section of the column as a^ 4 

ds 

dy 

■r-t and from the coordinates a. b of the still unknown 
as 

point 0. 



let us considor the roctangular flat plate of thick- 
ness 6, now assumed to bo constant, obtainod by cutting 
the column wall along a generating axis and developed in 
piano . 

If tho mass forces and tho external pressures acting 
on the plenes bounding the plate are everywhere zero, the 
unit stresses set up in the plate by any system of forces, 
acting in its median plane and however applied along its 
edges, may be expressed in the following form (reference l) 



v \ ( Bl z) = ^;cg - g> 

o z 

CT* B ( 8l z) = ^ 
Bs* 




(6) 



cr*a , cr* z are the normal stresses in direction of axes s 
and z, and T* is the tangential stress (axes z, s 
being disposed parallel to the edges of the plate of length 
2tt and I, respectively, and corresponding to the devel- 
opment of .the modian contour of one of the end sections of 
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the column and to the section separated along the generat- 
ing axis s = 0). 

If the condition is imposed which must be : • 



a* a (O.zJ = cr* 8 (2tt,b) 
(T* z (0,z) = ff z (2tt,1) 
T* (0,z) = T (2tt,z) 



(7) 



where "5 is^any value of z "between 0 and I, the 
functions a B (s,z), O" z (s,z), T (s,z) can be expanded 
in Fourier series. 

* — \ 

We con, for example, express a s (s,z) in the form 

of 

(j* z (s,z) = £ n An ("z) sin ns + Bjj (z) cos ns (8) 

where the values of the coefficients An ("5) e- n ^- Bn (z) 
depend on the value "z. Then the whole plate can be ex- 
pressed as: 

<T* a (s,z) = £ n - An (z) sin ns + B n (z) cos ns (9) 

with A n , B n as functions of z. 

Tho last expression twice intcgratod with respect to 
s givos: 

r(s,z) = s $ (z) + \|/(z) + £ n (- sin ns - ^ cos ns) 

> n n / 

(10) 

$ and v|/ both being functions of z only. 

Substituting this last expression in tho first and 
third of equation (6), wo obtain: 

+ " " /A" B" \ ^ 

o* B = s $ + \|/ + E n ( - ;;a a - sin ns - cos ns ) 

' Mil) 

/ A • B_ « \ I 

T* = + S n f - -J- cos ns + -J- sin nsj 
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The function a* a (o,z) must "bo poriodic with reapoct to 
a, whatovor the value of "5, whence equation (11) af- 

'f<j*a B : ~ ' ' * - - " ■ 

« (z) => H z.+ Z (12) 

trhoro H and k are oonetanta. 

function r(a,z) oDvioualy satiafloa (MaocwollU) 
dlfforontlal oquation: 



a* 1 + aljE + 2 r . = o 
a a* a z 4 aa a az« 



(13) 



ODtainahlo from the conaidoration of tho equilibrium of an 
element da dz of the plate. 

Substituting the valueB derived from equation (10 ) in 
thia last equation afforda: 

* I7 ( Z ) + 2n [(2An n - ^-J* - n 8 A^ ain na + 

IV 



^2B n » - - n B B n ) coa naj =0 (14) 



thia equation, being satisfied for any value of 8, should 
have at the aame time for any value of n (whole) and z: 



IV 



B 



IV 



- 2n a + n* ^ = 0 I 

- 2n a Ba" + n* ^ = 0^ 



(15) 



are : 



The general integrals of this differential equation 

V = M z 3 + H z a + P a + Q 
An = (oq + ^z) ooah nz + ( | n + "n n z) ainh nz >(16) 
Bn = (cp n + X n z) cd ah nz + (Y n + u) n z) ainh nz 
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whero the quantities M, S, P, Q, , P n , | n , T> n , (J^, Xn, 

Y n , u> n aro constants and depond upon the limiting condition 

which must satisfy the stress and strain condition of the 
plate. 



Such conditions are: 

o-J (s.O) = 0 
T* (s,0) = 0 

t* (s.l) = 0 
w* (s,l ) = - 



for z = 0 



for z=l 



} 



(17) 



(18) 



w(s) "being dofinod "by oquation (5), except that one con- 
stant which is unimportant need not ho determined. 

Prom tho first condition of equation (17) by virtuo 
of (9) and (16) follows: 



ota = 0 



«P n = 0 



} 



(19) 



and from tho second, while kooping in mind the second of 
equations (ll), (12), (16), and (9): 



H 


= 0 






in 


la 

n 




(20) 


*n 


= - *a 

n 







From tho Genoral relation: 

|1 . i fa . 5k) (M ) 

9b E \ a m / 

(E = Young's modulus, "J = Poisson's ratio) and from tho 
first condition of oquation (18) evaluated as hofore from 



the oipressions 0" 



and 



supplied hy equation (9) 



and from the first of oquation (11 ), we ohtain , while al- 



E.A..C..A. Technical Momoramdum ,Ha. -.8&4 



9 



lowing for (16), (9) and (20): 



*n 
*n 



- 3111 

*n ^n 
*n <" n 



(22) 



Hence , 



E n - 1 , m 

■ m n ■ 



2 + m + 1 i tiuah n i 
n m 



tanh n I - S-±-i I 
m 



(23) 



a quantity easily determined. 

Deducing equation (2) with respect to s, and then 
substituting the Becond complex derivative of t supplied 
by means of the derivation relative to z of equation (21), 
the formula assumes the general character of 



as 3 



1 9t 

o as 



1 ^£g 

e as 



_1_ 3g g 
In dz 



(24) 



Substituting in this last relation the values obtained 
for s = I, deriving the ratio of z (equation 9), the 
first part of equation (11 ), and the second part of equa- 
tion (11) with respect to s, and utilizing the conditions 
imposed by the Becond of- equation (18) the coefficients 
T) nt w n - and consequently, with equations (20) and (22) - 

the values of the coefficients Pn , Xn» £n> 8111 *n can De 
determined. 

Carrying out the operations indicated above and de- 
veloping equation (5) in Fourier series, putting 



dw(s) . 
ds 



- ^n (Jin 1 + l*n fl a + V*"*) sin n o + 



(l? ' + V" a + V n IM h) cos n a (25) 



we find: 
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M =0 



where 
*n 



J. D 

-*[( 



+ (v n ' + T? n " a + v n '" h) 
T) n _ - 

X = " + tV' a + tin'" h) 

a .. • ••• - - Tn ' — J 



3i + 3 x 5 



m 



m + 3 \ 
m n / 



sinh n I + 



+ ^3m + 3 



(26) 



U. I + £5a) cosh n i] (27) 
m n * J 



It is readily seen that the stress condition in the 
flat plate thus defined differs very little from the sec- 
ondary stresB condition in the column. 

The limiting conditions manifestly coincide, so that 
they are assigned to the state of stress of the plate rel- 
ative to the stress condition of the column. 

Voir consider a distance (or length) of the cylindri- 



cal wall of the 1 column of depth dz and length s 



o » 



denoting the angle of the tangent to the median contour of 
the straight section of the column at any point of the arc 
considered and the normal to it at the origin of this arc, 
the expression of tho equilibrium of displacement accord- 
ing to this last direction of the path of tho previously 
defined cylindrical wall, and assumedly constant thickness 
8 : 



(To 5 cos 



■ <■„> 



dz 



-/'if 



cos c (s ) ds dz (28 ) 



If, in this relation for 



and 



the values given 



in equation (6) for the flat plate are substituted, we 
will have after dividing both torms by dz : 



6 r— a c °s «(s 0 ) 



= 8 



J 



a 8 it 
asaz 8 



cos e(s) ds 



(29) 
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which, in reality, is identical with 




Ci). a f -f-£ 



a 



8 



CO 8 € 



TT CO 8 € ds-6 



T 3 a 3T d€ 

J as 8 ds 



sin e de 




r° a'T d£ r° O , . 

8 / r~ a — sin e da ■ / —7-^ sin € ds (31) 
J 3z s dS J p(s) 



o o 



whoro p(s) is the radius of curvature of tho cylindrical 
wall, tho approximation of equation (29) is so much great- 
er as the moan value of the ratio 8/p is less* 

In practice this ratio is generally quite small, hence 
the orror committed "by substituting for the secondary 
stress condition producod in the column wall as the result 
of rostrained wrinkling in one of the end sections, the 
corresponding ntress condition of tho flat plate, is prac- 
tically negligible. 

The stress condition of the column is given through 
equations (9), (ll), (16), (19), (20), (22), (23), (26), 
and (27) in function of the known quantity and of the 
etlll unknown coordinates a. b of 0. 

The last values can he determined by evaluating the 
condition which stipulates that the resultantt bending mo- 
ment of u 2 acting on the restrained section must be zero. 

Taking into consideration the relations (9), (16), 
(19), (20), (22), and (26), this condition can be expressed 
in the form: 



12 
a-n 

/ 
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x E n 7 n (I) [(n^i + i^" a + ji n "« b) cos ns - 

- ( v n J + V a + v n IM b) ain na] da 



> (32) 



y E n V n (I) Dl^n 1 + Hn" a + (in 111 "b) cos na - 

- (v n i + u n » a + u n MI "b) ain na] da f 

where 

E- I co ah n I + (l - ^a) a inh n I 

\ CD -— ^— ^ 



(33) 



From equation (32) we can now deduce the valueo a, b 
of the coordinates of 0, and ao determine the aecondary 
condition of atreaa of the column. 

Summing up, having obtained the values a, b from 
equation (32), and posting 



i n 



(34) 



= Vn' + V a + V" » 



the final expressions of the stresses existing in the col- 
umn can he written as follows: 
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0» = CT B = E - Ha 2 cosh nz +' Qt - ainh nzjcos na + 

+ v n * coBh nz + ^z - ainh nzj ain no 

CT B = a* = - + Rq a) coah nz + 

+ + z^ ainh nzj ain na + 1%^^ + a^ coah na+ 

+ ^ a + z) ainh nzj cob na 

T = ih + T * = ih " En Vn [ z COflh nB + 

+ ^H n z + ^ ainh nzj cob ns + Hjj. 
4* ^H n z + ^ ainh nzj ain na 



z coah nz + 

L 



(35) 



where the quantities , u n , and H n are given in equa- 
tions (34) and (23). 

So far, the thickness S has "been assumed constant 
hut by way of approximation the obtained findings may he 
oxtended to include the case of variable 6. 

Equations (6) and (13), which are the basis of the 
treatment, can be derived from the relations expressing 
the equilibrium of displacement according to b and to 
z of a plate element 8 ds dz, by substitution for the 
elementary stresses acting on each surface 8 ds, 8 dz 
of the elomont in question, the corresponding resultants, - 
and that is legitimate if the lattor lie in tho median % 
piano of tho plate. It is equivalent to considering o~ s , 
ct z , t 'as stresses applied to any linear element of the 
surface. 

However, if 8 varies so that the resultant of the 
stresses acting on each element 8 ds, 8 dz always lies 
in the median plane of tho plate, and if the effect of the 
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unit stresses acting perpendicularly to the latter can al- 
ways he kept negligible, the equations for the equilibrium 
of displacement of plate element ds, ds may still he 
written in the same form as for the plate of constant 
thickness, and equations (6) and (13/ can he made to retain 
thoir validity - it being, of courso, understood in that 
caso that o~* , a* z , and T aro stresses acting on olomont 

8 ds and 8 dz, respectively. 



EXAMPLE 



Consider a cylindrical tube having a straight section 
as in . figure 2; the unit length is chosen so that the perim- 
eter of the section is 2tt. The wall thickness 6 is con- 
stant. The length of the tube is I = 16. 

In the plane of the section a system of orthogonal car- 
tesian axes is established whorein axis x is the tangent 
to the ;.iodian contour of the section in one of the points 
where it intersects the axis of symmetry, axis y being 
coincldont with it (fig. 2). 

For evidont reasons of symmetry in points symmetrical 
with respect to y, the displacomonts and socondary stress- 
es assumo opposito values. It suffices therefore to ana- 
lyze only the bohavior of one-half of the soction; for ox- 
ample, that for which x > 0. 

Introducing in oquations (l) and (4), the values rela- 
tive to the present example, we have: 

M 

T = 0.19296 -=- 
6 

f| = 0.23395 

The section being symmetrical, it must be that a = 0; 

4 E supplied by equation (5), may be written: 
ds 

dw _ dwx ^ dw3 
ds ds ds 

dw x dw a 

The appended table I lists the values of ^jj— ^jj— 
along with those of the coordinates z, y of the median 
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dz dy 

contour and of ^ relative to the particular half of 

the . .aaction,- in function -of b. • 

In order to quickly determine the values of the coef- 

ficienta i^*, p^", n n l " , i? n » , v n ", i? n ,w , the method of 

Eunge and Emde (Of. G. Caesinis, Calcoli nuraericl grafici 
meccanici. Gap. XIV) is applied. 

Denoting with S«, S n ," D 1 , D" the sums and the dif- 

ferences of the values of — -, — - in points symmetrical 

ds ds 

with respect to y, tre have: 

1 w 1 * 

i V = 2 n S h » cos n s h - IV = D h « sin n s h 

| V" = lhS h » cos n s h i JV'"= |h Dh» sin n s n 

^h 1 = ^>." = ^ ^ or an y n reasons of symmetry, and 

since, i\o seen, a = 0 , we have: 

V = Hn 1 = Mta" = ^n"' = 0 
for any value of n. 

The execution of the calculation gives: 

V 0 ' •»£ = + 0.00233 W 0 "' = + 0.00088 

V -jf = + 0.17111 I?/ 1 ^ = - 0.27625 

V a f ^ = + 0.07340 U a ,M ^ = - 0.00878 

W B ' §L - - 0.02252 V 3 ,M fif = + 0.05433 

V 4 ' Y = " 0.03663 V4 fM ~ = + 0.00993...... 

u' = - 0.01159 U B ,M = - 0.01453 
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+ 



0.00564 



i it 



4f = + 0.00089 




+ 



0.00803 



11 



a + 0.00230 



+ 



0.00349 



s 



11 1 



= - 0.00292 



■ ■ ■ ■ t 



The values of the quantities R n , T n , and T n are de- 
termined "by means of equations (23), (27), and (33), put- 
ting n = 3.3 which , together with the first part of 
equation (32), give: 



Next, "by means of equation (35), we obtain the fol- 
lowing values of the stresses cr* z , T* for the restrained 
section (z = 16): (See table II.) 

T7o wish to thank Professor E. Fistolesi for the usoful 
advice ho has given during the working-up of this report. 



Translation by J. Vanier, 
National Advisory Cormittee 
for Aeronautics. 



b = 0.62031 
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TABLE I 



z . . . . 

y 

dz 

■ ■ • 

aa 
dy 

as - * ' 

as M 
dwa 06 
as y 



0 
. 0 

+1.00000 
0 

+.19296 
-.23395 



g 

+0.39270 
0 

+1.00000 
0 

+.19296 
-.23395 



TT 

4 

+0.78540 
0 

+1.00000 
0 

+ .19296 
-.23395 



3n_ 
8 

+1.13563 
+.15205 

+ .69284 
+.73C55 
+.02315 
-.15975 



2 

+1.26376 
+ .51173 

-.06743 
+ .99783 
-.11013 
+.01578 



5JT 
8 

+1.08812 
+ .85038 

-.77384 
+. 63192 
-.12186 
+.18104 



3TT 

4 

+0.75844 
+1.06376 

-.89776 
+.44046 
-.10862 

+.21003 



7TT 

_8 

+0.38929 
1+1.19522 

-.97411 
+.22608 
-.10001 
+.22789 



TT 

0 

+1.23980 
-1.00000 

0 

-.09709 
+ .23395 . 



TABLE II 



a 


0 


TT 

8 


TT 
4 


3TI 

8 


TT 
2 


5TT 

8 


3TT 
4 


7TT 

8 


TT 


* T* 5 
>V » ff* A 

u 


0 

-0.00803 


-0.01402 
-.00719 


-0.03087 
-.00264 


-0.02781 
+.00625 


+0.00801 
+.01254 


+0.02577 
+.00245 


+0.02347 
-.00264 


+0.01403 
-.00625 


•o 

• 

-.00765' 



These are the maximum values verified on the column; streBBes <s\ and T* disappear 
very rapidly from the point of fixation toward the free end. 
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Figs. 1,2 




Area = 2.5912 
Figure 2. 
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Jigs. 3,4 
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